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Abstract 

We propose a novel approach to incorporate the proton-proton (pp) Coulomb force into the three-nucleon 
(3N) Faddeev calculations. The main new ingredient is a 3-dimensional screened pp Coulomb t-matrix ob- 
tained by a numerical solution of the 3-dimensional Lippmann-Schwinger (LS) equation. We demonstrate 
numerically and provide analytical insight that the elastic proton-deuteron (pd) observables can be deter- 
mined directly from the resulting on shell 3N amplitude increasing the screening radius. The screening limit 
exists without the need of renormalisation not only for observables but for the elastic pd amplitude itself. 
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I. INTRODUCTION 



The inclusion of the Coulomb force into the analysis of nuclear reactions with more than 2 
nucleons is a long standing problem. The main reason is the long-range nature of the Coulomb 
force which prevents the application of the standard techniques developed for short-range interac- 
tions. One possible way to avoid the difficulties including the Coulomb force is to use a screened 
Coulomb interaction and to reach the pure Coulomb limit through application of a renormalisation 
procedure 

Mia. 

The problem is especially timely to be solved when considering the interaction of protons with 
deuterons below the pion production threshold. For this 3N system using the Faddeev scheme high- 
precision numerical predictions for different observables in elastic proton-deuteron (pd) scattering 
and for the deuteron breakup reaction are being obtained [5], however, only under the restriction to 
short-ranged nuclear interactions. The high quality of the available pd data for both processes re- 
quires, however, in the theoretical analysis the inclusion of the Coulomb force into the calculations. 
Furthermore the seminal progress [6f] in the development of nuclear forces in chiral effective field 
theory calls also for a very precise solution of the pd scattering equations to test unambiguously 
these new dynamical ingredients. This test can only be completely satisfactory if the pp Coulomb 
force is perfectly under control. 

For the elastic pd scattering first calculations, with modern nuclear forces and the Coulomb 
force included, have been achieved in a variational hyperspherical harmonic approach [?]]. Only 
recently the inclusion of the Coulomb force became possible also for the pd breakup reaction [§]. 
In contrary to [?J where the exact Coulomb force in coordinate representation has been used 
directly, a screened pp Coulomb force has been applied in momentum space and in a partial wave 
basis. In order to get the final predictions which can be compared to the data, the limit to the 
unscreened situation has been performed numerically applying a renormalization to the resulting 
on-shell amplitudes 0, Q]. This allowed for the first time to analyze high-precision pd breakup 
data and provided a significant improvement of data description in cases where the Coulomb force 
plays an important role [lo| . 

However, in spite of that substantial progress some important questions remained unanswered. 
One concerns directly the results of these calculations for two kinematically complete breakup 
geometries: the pp quasi-free-scattering (QFS) configuration, in which the not detected neutron 
is at rest in the laboratory system, and the space-star (SST) geometry, in which all 3 outgoing 
nucleons have the same momenta (magnitudes) in the plane which in the 3N cm. system is 
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perpendicular to the incoming nucleon momentum. The theoretical predictions based on nuclear 
forces only show, that the cross sections for QFS and SST are quite stable against changes of 
the underlying interactions, including also three- nucleon forces Q]. At energies below ~ 20 MeV 
theory underestimates the SST pd cross sections by « 10%, and overestimates the pp QFS cross 
sections by ~ 20%, respectively [a, |8|. With increasing energy the theoretical cross sections come 
close to the data, which indicates that the pp Coulomb force is very probably responsible for 
these low energy discrepancies. However, the Coulomb force effects found in 8[ are practically 
negligib 
data 
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e for the pc SST configuration and only slightly improve the description of the pp QFS 
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This inability to understand the pp QFS and pd SST cross sections calls for reconsidering the 
inclusion of the Coulomb force into momentum space Faddeev calculations. One main concern 
in such type of calculations is the application of a partial wave decomposition to the long-ranged 
Coulomb force. Even when screening is applied it seems reasonable to treat from the beginning 
the screened pp Coulomb t-matrix without partial wave decomposition because the required limit 
of vanishing screening leads necessarily to a drastic increase of the number of partial wave states 
involved. As an example we provide numbers for the exponential screening e~^R^". Taking the 
screening radius R = 20 fm and n = 4 requires all / < l max = 10 partial wave states to reproduce 
the 3-dimensional pp screened Coulomb t-matrix at E l p b = 13 MeV. Increasing the screening radius 
to R = 120 fm requires l ma x ~ 50 which is a big numerical challenge. Even more that would lead 
to an explosion of the number of 3N partial waves required for convergence. 

Another problem concerns the treatment of the pp Coulomb interaction in its proper coordinate, 
which is the relative proton-proton distance, throughout the calculations. Any deviation from this 
restriction can cause effects which are difficult to estimate. 

In the present paper we propose a novel approach taking both concerns into account to incor- 
porate the pp Coulomb force into the momentum space Faddeev calculations, in which we apply 
directly the 3-dimensional screened pp Coulomb t-matrix without relying on a partial wave de- 
composition. 3-dimensional solutions of the LS equation for different screening radii are used to 
approach the final predictions. We demonstrate, that even the physical elastic pd scattering am- 
plitude has a well defined screening limit and does not require renormalisation. In contrast in case 
of pp scattering the scattering amplitude requires renormalisation in the screening limit which, 
however, has not to be applied when only observables have to be addressed [if]]. 

In section [II] we present details of the formalism and in section IIIII the physical amplitude for 
elastic pd scattering. The screening limit is discussed in section [IV] and the results shown in section 
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PVT . The summary is given in section IVIl In various appendices detailed expressions for different 
kernel elements appearing in the Faddeev equation with the screened Coulomb force are included 
and different terms contributing to the on-shell elastic scattering amplitude are given. 



II. FADDEEV EQUATIONS WITH SCREENED PP COULOMB FORCE 

The 3N pd system can be regarded in a particle basis or in the isospin basis. The corresponding 
completeness relations are 

\pnp >< pnp\ + \ppn >< ppn\ + \npp >< npp\ = 1 , (1) 
|(*1/2)T - 1/2 >< (tl/2)T - 1/2| = 1 . (2) 

tT 

In ([2]) t is the 2-body isospin which together with the isospin ^ of the third particle is coupled 
to the total isospin T. We use the convention that the proton (neutron) has the isospin magnetic 
quantum number — \{\)- Since NN forces are usually given in an isospin basis we stick to that 
formulation also in the 3N system. Then a general NN force acts in the 3N system as 

V = ^2 Yl \( tl / 2 )™T >< (tl/2)TM T \V\{t'l/2)T'M T >< (t'l/2)T'M T \ . (3) 

t,T t',T> 

Here we assume that t is conserved, which in the pd system leads to [ijj 

V = |(01/2)l/2 - 1/2 > V™ < (01/2)1/2 - 1/2| 

+ |(ll/2)l/2 - 1/2 > (l/W™ + 2/Wp 1 - 1 ) < (11/2)1/2 - 1/2| 

+ |(ll/2)l/2 - 1/2 > ^V^ - Vj- 1 ) < (11/2)3/2 - 1/2] 

+ |(ll/2)3/2 - 1/2 > (2/3^4° + l/SV^ 1 ) < (11/2)3/2 - 1/2| 

+ |(ll/2)3/2 - 1/2 > ^{V™ - Vjjp 1 ) < (11/2)1/2 - 1/2| . (4) 
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We defined in the usual manner 



Vnp = < t0\V\t0 > , (5) 
VpV 1 = <1-W-1> • (6) 

The latter decomposes into the strong part and the pure Coulomb part, which we assume to be 
screened and parametrised by some parameter R 

yl-l _ ystrong , ycR fj\ 
pp pp ' PP ' v / 



Eq.(H]) exhibits that charge independence breaking leads necessarily to a coupling of T = 1/2 and 
T = 3/2 states and the complete action of the Coulomb force also requires the presence of both 
total isospin states. Exactly the same decomposition holds for the 2N t-operator t 



+ 
+ 
+ 
+ 



(01/2)1/2 - 1/2 > t™ < (01/2)1/2 - 1/2 



(11/2)1/2 - 1/2 > (1/3C + 2/34- 1 ) < (11/2)1/2 - 1/2| 
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I 1 I r> ) I /2 - 1 /2 > ^(t™ - t 1 ' 1 ) < (11/2)3/2 - 1/2| 



(11/2)3/2 - 1/2 > (2/3^° + l/St 1 - 1 ) < (11/2)3/2 - 1/2| 



1(11/2)3/2 - 1/2 > ^(t™ - t 1 - 1 ) < (11/2)1/2 - 1/2| . (8) 
It is not. difficult, to set' [17] tfiat tltis is consist. cut with the 2-body Lippmann Schwinger equation 



t = V + VG t 

when expanded into the 3N isospin basis. 

It is now convenient to define the three 3N isospin states 

| 7o > = |(01/2)l/2- 1/2 > 
I71 > = |(ll/2)l/2- 1/2 > 
| 72 > = |(ll/2)3/2- 1/2 > 



(9) 



Then t appears as 



(10) 



(11) 



7,7 



with 



tyrf — 8tt't t 



TT' 



(12) 



and tJ T ' can be read of from ([S]). 

We use the Faddeev equation in the form 



T|3> >= t>|$ > +tPG T|$ > 



(13) 



where P is defined in terms of transposition operators, P = P12P23 + -Pi3-P23; Go is the free 3N 
propagator, |$ > the initial state composed of a deuteron state and a momentum eigenstate of the 
proton. Knowing T|<1> > the breakup as well as the elastic pd scattering amplitudes can be gained 
by quadratures in the standard manner j^]. 



Expanding T|3> > 

T|$> = ^|7>^ 7 

7 

(fT3l) turns into 

r 7 = j^tyy < y|p|* > +X)*77' < vi-fGoiy > r 7 

In a more detailed notation we use (fT2"j) and define 

T 7 = if 

< 7 |P| 7 ' >= 5 TTI F{tt'T)P ttl = b T T<Ptt> 
where acts only in spin-momentum space, 

Ptt> = P12P23 + ( — ) t+t -P13-P23I spin-momentum 

and F(tt'T) is given by 



F(tt'T) = (-)* v / 2tTTv / 2t / + 1 



x 



t 1/2 1/2 

f r 1/2 



For the convenience of the reader the nonzero values are 



F(lll/2) = -\ 

F(113/2) = 1 

F(101/2) = ^ 

F(001/2) = ' 
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V3 



F(on/2) = 9 



This leads to the more explicit form of the coupled set of equations ([T5 

if = t f*p*i0 > + e tr£ ^ if . 

T' i" 

where we used the isospin independence of the free propagator Go and 



1$ >= |to > 10 > 



This set (|25|) would be the starting point for a calculation using momentum vectors instead of 
an angular momentum decomposition. In view of forthcoming calculations we provide that set as 
follows 



T i/2 _ V3,l 10 2 u , v3,l 10 2 f i-u p r T i/2 

~ 2 (3*^ + 3*pp 1 )- PllG!o - r i /2 + "3~(*np ~ *pp ^-PnCo^ 2 , (28) 



T i 3/2 - yg(*np ~ *pp 1 ) jP iol0 > +-^( t n P - *pp 1 )PioGqTI /2 

?(*np ~ ipp 1 )-?!!^^ 2 + (— + -tir 1 )PnG T^ / ^ . (29) 



3/2 n ^ pp ' 1 ^ 3 n P 3 



Note tp p 1 is driven by the strong and the screened pp Coulomb force ([7]) . 

Here, however, we apply the angular momentum decomposition and define the basis states 

\pqa >= \pq(ls)j(X^)I(jI)JM > (30) 

where p, q are the standard Jacobi momenta, (ls)j refers to the 2-body subsystem, (A|)J to the 
third particle and JM is the total angular momentum and its magnetic quantum number. In that 
partial wave basis the set (j25j) reads 

fx 11 
< pqa\Tj > = 2J / dp'p 2 dq'q' 2 < pqa\t t 2 \p'q'a' X p'q'a'\P t 2 \4> > 

a' 

+ E E / dp' P ' 2 dq'q' 2 < pqa\tJ T '\p' q 'a' X Mal\Fg,Gtf$ > . (31) 

T't" a' J 

This form allows us to separate the lower partial waves in the 2-body subsystem, where the strong 
force acts together with the Coulomb force, from the higher partial waves, where only the pure pp 
Coulomb t-matrix is present. 

Thus for the higher partial waves, denoted as \pqa > >, we obtain 

<pqa > \Tj > = (32) 
<pqa > \jf> = E / dpp 2 dqq 12 <pqa > \t T J\\p'q l a l ><pqa\P^\(t)> 

a' 

+ E E / dp'p ,2 dq'q' 2 < P qa>\tJT'\p'q'a' X p' q' a'\P^G Q \Tj,' > (33) 

T't" a' ^ 



where tf c T ' up to a factor is given by the screened pure Coulomb t-matrix contained in tp p 1 . 
Thus 

t TT = fvnrg (34) 

and f 1T1T ' can be read of from © as f 1 ^ = § = 2/ 1 f 1 f,/ 1 5 1 f = J 1 ! 1 ! = . 
Inserting the factors f 1T1T it simply follows from (|33|) that 

< pga > |T ] 5 > + v / 2 < pga^T^ >= . (35) 

The physical meaning of that relation is trivially evident in the particle picture. One has 

[*yi > = —=(\ppn > +\pnp > —2\npp >) (36) 
V6 

1 72 > = —j={\ppn > +\pnp > +\npp >) (37) 
v3 

where the sequence pnp for instance corresponds to the sequence of particle numbers 1, 2, 3. 
Inserting the definitions (|14p and (|16|) the relation (|35|) can be rewritten as 

<pqa > \(<^i\ + ^2 <y 2 \)T >= J^<pqa > \(<ppn\+ <pnp\)\T > . (38) 

Now the T-amplitude projected on high angular momenta is proportional to the Coulomb pp 

t-matrix, which acts in the 23-subsystem. Consequently the overlap in f|38|) is zero. 

i 

Now using (|35p and (|34p the relation (|33[) for < pga > |T 1 2 > simplifies to 

< pqa>\T^ > = ~^J2 [ dp'p' 2 dq'q' 2 < » > '\tpp\p V 'a >< p q a \Pi \<j) > 
V3 , J 



+ 4=E / d PP' 2 dq'q' 2 < pq^ltpplp'q'a >< p'q'a\P w G \T 2 > 

a' 

- Je/ *VW 2 E / d*riPtf^<1>q* > \*gW > 

a' a" 

< p'q'a!\P n G () \ P "q"a' < >< tf'tl'a^T? > +v^|T^ >) . (39) 

Because of (|35|) there are no contributions of the higher partial wave components and only the 
lower ones, denoted as < pqa < \, are applied onto the amplitudes \T >. Due to ([35]) the same 
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is true for < pqa > \T-^ >. For the sake of a simpler notation we introduce \A >= \pqa < > and 
\B >= |pq , a > > in the following and drop the summation and integration signs. 

Now we turn to the lower partial waves. From (j27j) - (f29j) and inserting , P^, explicitly we 

get 

<A|T 1/2 > = -l<^| t oop oo |0> 



> 
> 



-\< A\t™P 00 G \A' >< A'\Tl /2 > 

- ^ < A\t™P 01 G \A' >< A'\Tl 12 > 

- ^ < A\t™P 01 G \B' >< B'\t} /2 > , (40) 
<A\Tl /2 > = -L < A\(t™+2t 1 ^ 1 )P 10 \ct>> 

+ -L < + 24; 1 )P 10 G |^ X A'|T 1/2 > 

- i < A\(t™ + 2^ 1 )P 11 G |A / >< A'\Tl /2 : 

- i < A|(t£ + 24- 1 )P 11 G | J B / X P'lT^ 2 
+ ^ < - ^AiGblA' >< A'|if 2 > 
+ ^ < A\(t% - t^PnGolB' x P'llf 2 > , (41) 

< A\T* /2 > = -L < A|(4° - tJ^Piol* > 

+ -±= < A\(t™ - ^PiqGqIA' >< A'\Tl' 2 > 

- < - t^^PnGolA' >< ^'IT, 1 / 2 > 
+ i < A|(2t£° + t^PnGolA' >< A'\T* /2 > 

- < A|(t£ - ^PnGolP' >< P'lT^ 2 > 

+ < A|(|<£ + i^-^PnGolP' X P'lTf/ 2 > . (42) 

Note tpp 1 is the pp t-matrix driven by the strong and screened pp Coulomb force and is projected 
on a certain set of low partial waves. 

In (|4"T|) and ([15]) one can use ([55]) to eliminate < P'|T 1 3/ ' 2 > leading to 

<A\Tl /2 > = -L<A\(fij, + 2t 1 - 1 )P 10 \<l>> 

+ -L < A|(t£ + 2^- 1 )P 10 Go|^ / X A'|T 1/2 > 

- 1 < + 2ti; 1 )P 11 G |^ / X A'^ 2 > 
+ ^ < - t^PuGolA' >< yl'lTf/ 2 > 

- 1 < A|^PnGo|B' >< 5'IT^ 2 > , (43) 

< Airf 72 > = -L < a|(4° - t^ 1 )^ > 

+ < A|(t™ - ^^PwGolA' x A'\T^ /2 > 

- -L < A\(t 1 n p -tl- 1 )P 11 Go\A'><A'\T^ 



nl/2 > 
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+ - < A|(2t£° + t l - l )P n G \A! >< A'\T? /2 > 

- 4= < ^l^pAiGol^' >< S'lT^ 2 > . (44) 
V2 

The sum over the high partial waves acts on |T X > which, according to (|39p . is driven by the 
screened Coulomb t-matrix t%£. In order to avoid approximations we sum up the high partial waves 
to infinity and add and subtract the projection of tpp on the finite number of low partial waves. 
Thus we put 

\B'>< B'\T{ 12 >= \fl /2 > (45) 



where according to ([39 



\fl /2 > = -^JPi(# > +-^t c *P 10 G \A' X A'\T* > 

- \t£PuG \A' X A'(\T? > +V2\t} >) (46) 



with 



C = C - \ A >< a Kp\ a ' >< A '\ • ( 47 ) 

In (I46p the projection on partial waves to the right of tp 1 ^ has been carried out to infinite order 
leading to the 3-dimensional screened Coulomb t-matrix tpp. But due to (|47ft also the partial wave 
projected matrix occurs. Now we insert (14"6|) into the set (0Q|), (03]) and 

<a\tI' 2 > = ~<A\t^P^\4» 



i < A\t™PooGo\A' >< A'\T^ > 
/3 



v 3 < A\t™P 01 Go\A' ><A'\Tl /2 > 

<A\t n °pP 01 Go\fl /2 > , (48) 



^3 | o n../n,.i^i/2 
2 ^ -^rnp 



<A|T 1 1/2 > = -±= < A\(t 1 n °p + 2t 1 pp 1 )P w \ ( f ) > 

+ -L < A|(t£° + 2tir 1 )P 10 G |A / >< A'|T 1/2 > 



2^3 IVnp pp 

i < A\(t™ + 2t p ; 1 )P 11 Go|^ >< ^'ir^ 2 > 

2 /1 H t 'np 6 pp 



1^ 4|,10 r,.,< 
2 



- - < AI^PhGoIT^ 2 > , (49) 



lllf 2 > = -L < - t^Pio^ > 
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+ ^=< A\{tnp ~ &PyoG Q \A X A'\T l J 2 > 



+ i < A|(2t£ + ^^PiiGol A' >< A'llf 2 > 
1 

3 



< AI^PnGoir^ 2 > . (50) 



Eqs. (|48p - ([50]) together with (j36]) form a closed set for amplitudes projected on low partial 
waves only. Since we anyhow iterate the equations and sum up the resulting series by Pade one 
can stay with that form. However, we prefer to eliminate (|46|) to arrive at our final form of coupled 
equations: 



<A|T 1/2 > = -\<A\t™P m \$> 
1 



2 



- < A\t^ p P iG tpp P 10 \(j) > 



2 

\ < A\t™PooGo\A' >< A'|T 1/2 > 

\ < A\t™P Q1 G i c *P 1Q G G \A' >< A'\t} > 



^ < A\C p P 01 G \A' x A'\Tl /2 > 
1 



2 



+ ^= < A\t™P Q1 G t c Xp n G \A' x A'\T* > 
+ ±= < A\t™P 01 G t c Xp u G \A' x A'\T* > , (51) 



<A\T l J 2 > = -^ < A\(t 1 n ° p + 2t 1 p ; 1 )P w \^> 
1 



2^3 



2^ 



+ -L < A|(t£ + 2tJ- 1 )P 10 G |A x A'|T 1/2 > 

- ^= < AI^PhGo^PxoGqIA x A'|t| > 

- i < A|(t™ + 2ti; 1 )P 11 G |A x A'\tI /2 > 
+ I < AltMPnGotgPnGolA' >< A'\T* > 
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+ ^ < A|(45 - ^^PnGolA' x A'\T^ 2 > 



+ ^ < A\t™PuG t$PnGo\A' >< A'\T* > , (52) 



<A\Tl ,2 > = ±=<A\{t%-t^ l )P 1Q \j>> 
- —f= < A|t^pPnG o tppPio|0 > 
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+ -±= < A|(t*° - ffiPnGolA' >< A\T l J 2 > 
-4 < A\t™P n G t c Xp 10 G \A' >< A'|T ' > 



V6 



1 < - tl-^PnGolA >< A\Tl' 2 > 



3^/2 

+ i < ^|(2ti° p + ^ 1 )P 11 G |A' >< A'|2? /2 > 

+ i < AI^PuCoi^PnGolA' >< A']!? > . (53) 

Note that in all three equations there occur terms where the full 3-dimensional screened Coulomb 
t-matrix is sandwiched between permutation operators and occurs together with the strong np 
t-matrix. Due to (47) there are also other terms of that second order type in t-operators where, 
however, is projected onto the lower partial waves states. 

Whereas all the other expressions are of our standard type [5[ the ones with the 3-dimensional 
tpp operator require special care and will be dealt with in appendices. 

Now we would like to write down equations (|5ip , (|52ll and (|53]) directly in our standard momen- 
tum space partial wave basis \pqa > which is an extension of \pqa > by adding isospin quantum 
numbers 

\pqa >= \pqa > |(t~)T >= \pq(ls)j(\±)I(jI)J(t±)T > . (54) 

Though this appears as a repetition it is a necessity to cast the set into our standard form, which 
underlies our existing codes. 

We distinguish between the partial wave states \pqa > with total 2N angular momentum j 
below some value j max : j < jmax, in which the nuclear, Vn, as well as the pp screened Coulomb 
interaction, Vj 1 (in isospin t = 1 states only), are acting, and the states \pqf3 > with j > j max , 
for which only Vj 1 is acting in the pp subsystem. The states \pqa > and \pq/3 > form a complete 
system of states 

/ p 2 dpq 2 dq(^2 \PQ a ) (PQ a \ + \ pq ® = L ( 55 ) 

ot f3 

Projecting Eq. (|13|) for T|$ > on the \pqa > and \pqj3 > states one gets the following system of 
coupled integral equations 

(pqa\T\$) = (pqa\t% +c P\<5>) 



+ (pqa\t% +c PG J2 J p' 2 dp'q' 2 dq'\p'q'a')(p'q'a'\T\<5>) 
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+ {pqa\ t* +c PG / P ' 2 dp'q ,2 dq' \p'q'(3') (p'q'^T |<&> 



(56) 



(pqP\T\$) 




(57) 



where t^ +c and are t-matrices generated by the interactions V/v + Vj 1 and V c r, respectively. 

The third term on the right hand side of (|57p is proportional to < pqf3\t^PGo\p'q' f3' >< 
p'q'j3'\t^. A direct calculation of its isospin part shows that independently from the value of 
the total isospin T it vanishes. This corresponds to the result found in (I35p . 

Inserting < pq(3\T\& > from (|57l) into (f56l) one gets 



We used again a relation corresponding to now in isospin notation. This is a coupled set 

of integral equations in the space of only the states \a >, which incorporates the contributions 
of the pp Coulomb interaction from all partial wave states up to infinity. Using the definition of 
the 1 7 >i states from (fT0|) a direct calculation shows that the set (f58|) is identical to the set (f5Tj) . 
()52p and ()53p . It can be solved by iteration and Pade summation. There are two new leading 
terms < pqa\t^ +c PG t^P\^ > and -< pqa\t^ +c PG \a' >< a'\t^P\^ >. The first term must 
be calculated using directly the 3-dimensional screened Coulomb t-matrix t^, while the second 
term requires partial wave projected screened Coulomb t-matrix elements in the \a > channels 
only. The kernel also contains two new terms. The term < pqa\t^ +c PGot^ PGo\a' >< a'|T|$ > 
must again be calculated with a 3-dimensional screened Coulomb t-matrix while the second one, 
-< pqa\tpf +c PGo\a' >< a'\t^PGo\a" >< a"\T\& >, involves only the partial wave projected 
screened Coulomb t-matrix elements in the \a > channels. The calculation of the new terms with 
the partial wave projected Coulomb t-matrices follows our standard procedure. Namely the two 



(pqa\ r|#) = (pqa\ t% +c P\<S>) + (pqa\ t« +c PG t«P\$>) 

~ (pqa\4 +c PGoY^ I p 2 dp'q' 2 dq'\p'q'a')(p'q'a'\tKp\<S>) 




(58) 
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sub kernels t^ +c PGo and t^PGo are applied consecutively on the corresponding state. The details 
how to calculate the new terms with the 3-dimensional screened Coulomb t-matrix are given in 
Appendix |A"1 

Now we turn to the physical scattering amplitudes. 

III. THE ELASTIC PD ON-SHELL TRANSITION AMPLITUDE 

The transition amplitude for elastic scattering is given by [j], [3] 

($'\U\$) = (^'IPGq 1 +PT\$) . (59) 

The amplitude for elastic scattering given in (|59p has two contributions. The first one is independent 
of the pp Coulomb force 

< Q'IPGq 1 ^ >=< &\PV\$ >= -- < (/>'\PooV™\<f> > . (60) 

For the second one we use (|17|) . (23) and (24) and decompose into lower and higher partial waves 
which yields using (46) 

< t>>\PT\t> > = -\ < <mo|T 1/2 > < <P'\Poi\Tl /2 > 

= ~< 4>'\Poo\A X A\T l J 2 >-^-< <t>'\P Q i\A X A\T l J 2 >-^-< <j>'\P Q i\fl /2 > 
= -\ < <f>'\Poit c *Pio\<t> >~\< <j>'\Poo\A >< A\Tl /2 > 



X - < ^Pbxt^PioGolA >< A\Tl 12 > ~Y ' ^ p ~ 1 A ^ Mr}/2 



+ ^= < <p'\P 01 t^P u G \A >< A\{\Tl' 2 > +v / 2|T 1 3/2 >) . (61) 

Again besides standard low partial wave contributions terms occur with the 3-dimensional screened 
Coulomb t-matrix sandwiched between permutation operators. 

To calculate the matrix element (|59p one needs (pq\T\<&) composed of low and high partial 
wave contributions for T|<& > (see (57)). Using the completeness relation f|55[) one gets: 

(pq | T |$) = {pq | ^2 / p' 2 dp' q' 2 dq' \p'qa) (j)'q'a'\T |<J>) 
a' ^ 

~ (PQ I ^2 J p' 2 dp'q' 2 dq \p'q'a) (jp q a P |<&) 

a 1 

J P' 2 dp'q' 2 dq' \p'q'a')( P 'q'a'\t*PGoY^ / p"Vq" 2 dq" \p"q"a") (p"q"a"\T |$) 

a' a" 

+ (PQ I t*P\$) + (pq I tfPG J p' 2 dp'q' 2 dq' \p'q'a') (p' q' a' \T \<S>) . (62) 

a' 
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It follows, that in addition to the amplitudes < pqa\T\<& > also the partial wave projected 
amplitudes < pqa\t^P\& > and < pqa\t^PGo\a' >< a'\T\$> > are required. The expressions 
for the contributions of these three terms to the transition amplitude for the elastic scattering 
(and the breakup reaction) are given in Appendix [Bj The last two terms in (|62p must again 
be calculated using directly the 3-dimensional screened Coulomb t-matrices. In Appendix O the 
expressions for (pq\t^P\$>) (breakup) and (pq \ Pt^P |<&) (elastic scattering) are given together 
with the expression for < Q'\PGq^~\Q >. In Appendix [D] we describe how to get the last matrix 
element < pq\t^PG \a' X a'\T\® >. 

IV. THE SCREENING LIMIT 

The set of coupled Faddeev equations (|58p or (|48 |) -(|5l) |) is well defined for a finite screening radius. 
It is an exact set assuming that the strong NN t-matrix can be neglected beyond a certain j max , 
which is justified. Further the pp screened Coulomb force is taken into account to infinite order 
in the partial wave decomposition in form of the 3-dimensional screened Coulomb t-matrix . 
The important challenge is to control the screening limit for the physical on shell amplitude (|59p . 
The contribution (|60p is well defined and independent of the Coulomb force. The corresponding 
expression without partial- wave expansion is given in (|C5p , 

In the second part given in (|6ip there are contributions without and with an explicit t c ^ operator. 
The first term in (|6ip is responsible for the Rutherford scattering between proton and deuteron. 
It has the form 



< 0'|P O i&io|0 > = < 0'|PoiC p iol0 > - < ^01 \A >< A\t c *\A' >< A'\P 10 \<P > . (63) 



The 3-dimensional screened Coulomb t-matrix enters that matrix element totally off-shell, except 
for forward scattering. To show that it is sufficient to regard only the momentum space part of that 
matrix element. Using well known relations among Jacobi momenta and the symmetry property 
of the deuteron one obtains 



With Pqi = P w = P12P23 — -P13-P23 and using simple symmetry properties one obtains 



< 0'|PoiC p iol^ > 



2 < cj)'\P 12 P 23 t c ^(P 12 P 23 - P 13 P 23 )\<P > 



(64) 
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with 

It is easy to see that (q ' + l/2q') 2 /m ± E(q') ± (q + l/2q') 2 /m and obviously (q ' + l/2q') 2 /m / 
(qo + l/2q ') 2 /m except for forward scattering. Now the off-shell screened Coulomb t-matrix 
converges for R — > oo towards the pure off-shell Coulomb t-matrix, which is well defined. Thus 
except for forward scattering that matrix element has a well defined screening limit. 

It is interesting to regard the replacement — ► V£*f = ^re~ r l R which in the limit R — > oo 
turns ([65]) into 

-e 2 r 

<<t>'\P m V$P lQ \<t>> ~ — / d 3 q'<P d (l/2qo ' + q')Ml/2q Q + q') 

{qo - qo Y J 

- I d'q'Ml/Zqo ' + q')MU^0 + q')——^—^ ■ (67) 

J {qo + qo + q Y 

The first part is the Rutherford scattering amplitude multiplied by a deuteron form factor and the 
second term a correction due to antisymmetrisation of the two protons. 

The partial wave projected piece in (j63|) has also a well defined screening limit due to the same 
reasons. 

Next we regard the two matrix elements of in (61) integrated together with \A >-projected 
T-amplitudes. As before is off-shell from the left. On the right there is an intermediate 
momentum integration, which includes the on-shell point. In that case t£r acquires an infinitely 
oscillating phase factor in the screening limit [19] which, however, is integrable. Thus it remains to 
verify that the free propagator singularity of Go does not coincide with that on-shell singularity of 
tpp . Inserting the momentum representation and working out the permutations the free propagator 
singularity is located at 

£ 3 3 2 = 

4m m 4m 

where q ' and q are integration variables. On the other hand the singularity of t^{p,p 1 , E{q')) 
arising for p — p ' = is located in the matrix element either at qo ' + q ' + q = or at % 1 = q. 
This inserted into (|68|) yields in both cases a non vanishing expression. Thus in that respect all 
the three terms in (|6ip containing have a well defined screening limit. It remains to consider 
the partial wave projected T-amplitudes < A\T^ >, which enter in the matrix elements in (|6ip . 
They are the solutions of the coupled set of Faddeev equations (|58|) which are ill defined in the 
screening limit. However, we expect that in that limit these amplitudes < A\Tf" > acquire only 
infinitely oscillating logarithmic phase factors, which can well be integrated over. In fact as will be 
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demonstrated in a forthcoming paper only the on-shell amplitude < A\T^ > with p 2 + ^q 2 = mE, 
which appears in the breakup transition amplitude, acquires that oscillating factor. 

Thus we come to the conjecture that the physical on-shell elastic pd amplitude has a well defined 
screening limit and does not require renormalisation. This might appear surprising at a first glance, 
since the on-shell pp scattering amplitude has not a well defined screening limit. Our explanation 
is that in the pd system the Coulomb force does not act between a proton and the center of mass of 
the deuteron but between the two protons, where one of them is part of the deuteron and therefore 
the pp Coulomb t-matrix is integrated over the deuteron wave function in the final state in the 
elastic scattering amplitude. 

That analytical insight is well supported by our numerical results laid out in the following 
section. 

The case of the pd breakup process is quite different and will be dealt with in a forthcoming 
study. 

V. NUMERICAL RESULTS 

To demonstrate the feasibility of our approach we applied the outlined formalism to a simple 
dynamical model in which the nucleon-nucleon force was restricted to act in 1 S'o and S S± — 3 D\ 
partial waves only and taken from the CD Bonn potential [20]. The pp Coulomb force was screened 
exponentially 

V p R p (r) = (69) 

with the screening radius R and n = 1. 

To investigate the screening limit R — > oo we generated a set of partial-wave decomposed t- 
matrices, t^, based on the screened pp Coulomb force alone or for t^ +c combined with the pp 
nuclear interaction, taking R = 20, 40, 60, 80, 100, 120 and 140 fm. With that dynamical input 
we solved the set of Faddeev equations (|58p for the total angular momenta of the p-p-n system 
up to J < and both parities. In this first study we restricted ourselves to the perturbative 
approximation for the 3-dimensional pure Coulomb t-matrix: = V c , . Of course in the future 
studies that approximation will be avoided and the full solution of the 3-dimensional LS equation 
for the screened pp Coulomb t-matrix will be used. 

In Fig. [1] we show the convergence in the screening radius R of the pd elastic scattering cross 
section and compare the pd and nd elastic scattering angular distributions at the incoming nucleon 
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energy = 13 MeV. On the scale of the figure the cross sections for R = 40 — 140 fm are 
practically indistinguishable. The detailed picture of that convergence is depicted in Fig. [21 where 
the ratio of the cross sections obtained with the screening radius R to those with R = 140 fm is 
shown as a function of the cm. scattering angle G c . m .. It is clearly seen that already with the 
screening radius R = 40 fm converged results for the cross section are achieved. Increasing further 
the value of R provides cross sections which differ less than ~ 1% up to the forward scattering 
angles Q c .m. ~ 10°. At very forward angles, where the pp Coulomb force is dominant, larger 
screening radii are required to get the cross section with the same precision. 

The angular distributions shown in Figs. [T] and [2] were obtained taking in the elastic scattering 
transition amplitude (|59p the exchange term ($>'\PGq together with the first four terms in (|62p 
contributing to ($'|PT|$). In Fig. [3] we present how each term contributes to the cross section. 
When all terms are taken into account the resulting angular distribution is given by the solid line. 
The ($'|PT|$) term ( dotted line related to the first term in (|62p ) contributes significantly at all 
angles. At backward angles the largest contribution comes from the exchange term ($>'\PGq 
(dashed line) while at forward angles the most important is the "Rutherford" term ($'\Pt^"P\$) 
(double-dashed-dotted line related to the fourth term in (|62p ) calculated with the 3-dimensional 
screened Coulomb t-matrix tf' (in this first study treated perturbatively as tf' = V^). The two 
terms based on the partial-wave projected Coulomb t-matrix, (§'\PtfP\<&^ (dashed-double-dotted 
line related to the second term in (f62j) ) and {^'\Pt^ PGqT\<&^ (dashed-dotted line related to the 
third term in (|62p ), are about 2-orders of magnitude smaller and thus of minor importance. The 
fact that at very forward angles the contribution of the (<F '\Ptf PGqT\^ is an order of magnitude 
smaller than the contribution of the ^'\Pt^P\^ seems to justify the neglection of the last term 
(&\Ptf}PGoT\&^ in (|62p coming with the 3-dimensional screened Coulomb t-matrix. In future 
studies this term will be calculated to verify this statement. 

In Fig. 2] we demonstrate numerically that the elastic pd amplitude has a well defined screening 
limit and does not require renormalization. The real and imaginary parts of the partial wave 
contribution (&\P(T — t^P — t^PGoT)\^ to the elastic transition amplitude are shown for two 
combinations of the incoming and outgoing deuteron and proton spin projections and a number 
of screening radia R = 20, 40, 60, 80, 100, 120, and 140 fm. The additional term (|60p is real and 
independent of the screening radius. The fourth term in (|62p is also real under our approximation 

= and for angles different from zero has a well defined screening limit. Moreover it is peaked 
in forward direction and would dominate terms shown. All lines are practically overlapping. That 
shows that not only the cross section but the pd elastic scattering amplitude itself does not develop 
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an oscillating phase in the infinite screening limit. 
VI. SUMMARY 

We developed and presented a novel approach to include the pp Coulomb force into the momen- 
tum space 3N Faddeev calculations. It is based on a standard formulation for short range forces 
and relies on a screening of the long-range Coulomb interaction. In order to avoid all uncertainties 
connected with an application of the partial wave expansion, unsuitable when working with long- 
range forces, we apply directly the 3-dimensional pp screened Coulomb t-matrix. Furthermore we 
strictly insisted to treat the pp Coulomb force in its proper coordinate. 

Using a simple dynamical model for the nuclear part of the interaction we demonstrated the 
feasibility of that approach. We provided analytical arguments and showed numerically that the 
physical elastic pd scattering amplitude has a well defined screening limit and therefore does not 
require renormalisation. Well converged elastic pd cross sections have been achieved at finite 
screening radii. In this first study the 3-dimensional screened pure Coulomb t-matrix was replaced 
by the screened Coulomb potential and only a small number of partial wave states for the NN 
interaction was taken into account. This restriction will be removed in a forthcoming article. 
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APPENDIX A: MATRIX ELEMENTS WITH 3-DIMENSIONAL SCREENED 
COULOMB T-MATRIX: < pqa\t^ +c PG Q tf PG a \a' >< a'|T|$ > AND < pqa\t% +c PGot*P\§ > 

These terms must be calculated using directly the 3-dimensional screened Coulomb t-matrix t^. 
Introducing partial wave states \pqf3 > in LS-coupling 



\pq(5) = 



pq{l\)L{s\)S{LS)J(t\)T 



and using the recoupling 

\pqa) = Y, 



LS 



I S j 

A \ I 
L S J 



the matrix element (pqa\ t^ +c PGot^PGo \a') (a'\ T |$) can be written as 

{pqa\ 4 +c PG t?PG \a') (a'\T I*) = E / ^ C^fe^ 9)^0 (p, E(q)) 

a 

Y (&\P) Y (^V) Yl ( L ^ SM - /"I JM ) (L'v'S'M - n'\JM) 

a', 13' n,n' 

j p' 2 dp'q' 2 dq'G»(p\q')Z(pqp'q') 



where E{q) = E — ^-q 2 , G (p,q) = p2 1 , : and 



m 4m 



Z(pqp'q') = S S S'S,j l f Jl 'S MT _ 1 / 2 S MTl _ 1/ 2^ytt 



- J 1/2 1/2 4 
1/2 T 1 



> < 



1/2 1/2 1 
1/2 T f 



(-1) 



l+s+t' 



ss' < 



1/2 1/2 s 
1/2 5 s' 



X 1 (pqP,p'q'P') + 5 ss/ X 2 (pqP,p'q'P') 



X\ and A2 are given by 



(Al) 



(A2) 



(A3) 



(A4) 



X 1 (pq(3,p'q'(3') = (pq(l\)L» P™P™ 3 t?(T,T>)P™P™ \p>q> (l> \>)L> M '> , (A5) 



X 2 (pqP,p'q'(3') = (pq(l\)Lp J P%P™t?(T,T')P£P%\p' V '(I'X'W >') , 
and the isospin matrix element of the 3-dimensional screened Coulomb t-matrix is 



(t-)Tm T 



(t'-)T'm T >) = S tl 5 w S MT M T ,S MTir _itf(T,T') 



(A6) 



(A7) 
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with 

${T= l ~,T'= l -) = 2 -t* (A8) 

£(T = \>? = \) = t*{T = ^T' = \) = (A9) 

t?(T= 3 ~,T>= 3 -)= 1 ~t? (AID) 

The transposition operator P^ has been decomposed into three parts: acting on isospin-, P?- 
on spin-, and P™ on momentum-components of the basis |pga > 

p _ nm pm ps ps p* p* i pin pm p« ps pi pi / A 11 \ 

^ — ^12- r 23- r 12- r 23- r 12- r 23 ^13 ^23 A 3 *23 -<13 *23 • l Aii J 

Because Xj are scalars diagonal in L and // and in addition independent on [i one can write 
X x = ^p q (l\)L t M\pSP^t^(T : T , )P^\p , q'(l , \')L' tJ ,') 

= 4l7V\ ^ ( Pq $ X)LfJ \ P i m 2 P ^c(T,T')Pr 2 PS \p'q'(l'\')Lv) . (A12) 

Inserting in (|Al2j) from the left and right a complete system of states \pq > one gets 

Xi = J dpipl J dpi 'p\ ' J dqiql J dp J dq J dp ' J dq ' J dpi J dpi J dqi 
^^^^)tf(T,T';p 1)Pl >,p 1 op 1 ';E(qi))Y^(p \q') 

8(0- Qi) s (Pi + Q+ 2^i)KPi 2?1 ~ q')^(p' + ?i + ') • ( Al3 ) 

The quantity Yl Y~^(p, q)Y^,{p ', q ') is a scalar and therefore can only depend on scalars formed 
from p, g, p and g '. Only 5 of the possible scalars are independent e.q. p o g, p o p p o pq 
q o p ' : and qo q and they determine the geometry of p, q, p ' and g '. Using the 5- functions in 
(|A"l3l) it follows that 

1„ ^ 
P = 2 g + ?1 
1 _ 

Pi = -5 - 2^1 

-> / 1- 

Pi = g-ffi + 9 

p' = -&-\q' • (A14) 

The above 5 scalar products can be expressed by the only independent scalar pi o p 1 ' and by the 
magnitudes of seven vectors: p, q, p <f ', pi, q\ and pi '. This together with (|A14|) enables to 
perform the integrations in (|A13P by choosing the vector qi\\z and q £ (x — z). Such choice of 
coordinates makes that also p and pi £ (x — z) and: 4>p = 0, 4> Pl = vr, <p p i = (p q i, <jfy = <p q i + it. 
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As a result of the integration one gets 

-max 

2 



pp'qq' 



2tT 

| d<p q ,t?(T,T'; yp 2 + jjg 2 - jjg 2 , /2 + ^ /2 - jjg 2 ,Pi °Pi ',£(gi)) 



with pi o p' x = cos9 Pl cosOp'^ + sinO Pl sin9 p > cos(f> q ' and 



COS0. 



Pi 



cosO pi i 



<ii\l f> 2 + h 2 - hi 

T) /2 - In 2 _ I« /2 
P 2^1 4^ 



In a similar way one gets for X2 



(A15) 



(A16) 



X 2 (pqf3,p'q(3') 



32tt 2 
pp'qq' 



dq 



(-l)^LL' 



2L + 



27T 

| d^(T,T'; sjp> + ^ 2 - ^/p ' 2 + \q' 2 - \ql -pi oft ',£7( 9l )) 



Y V\i{6p 'Av' =<t>q' +'*,0q'Aq') ■ 



(A17) 



In (|A15j) and (IA17j) the limits of integrations over gi are 



g^ 111 = max 



P ~2 



P ~2 



(Ai8) 



The limits (|A18j) restrict the region of p'-values at given p, q and q' , for which Xi and X2 are 
nonzero. For g' < 2\p— 1| this requires that p' € (|p — f | — f-,P + § + > for 2|p — 1| < q' < 2(p + 
|): p' E (o,p+f + ^Y and for g' > 2|p- || and g' > 2(p + §): p' E - (p+f),£ +P+| 
Inserting Xi and X2 from ()A15p and (|A17|) into (|A3|) and (|A4p and using the orthogonality of the 
CG-coefhcients one gets finally 

{pqa\4 +c PGot?PG \a')( a '\T\$)=Y / [ P 2 dp t^ c ' R (p,p; E(q))G (p,q) 

-j 1/2 1/2 t 1 f 1/2 1/2 1 
1/2 T 1 j 1 1/2 T t' 



^2 (a (P' I Q ') S SS'6w'tiMT-i/2 S M T ,,-l/26Vtt 



a', 0i 
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/ 



p' 2 dp'q' 2 dq' 



(-i) 1+s+t '^\ 1/2 1/2 s \x 1 + s ss/ x 2 

1/2 S s' 



G (p',q) <pqa\T\<5> > (A19) 



Performing the analogous steps and using 



(p'qV | $} = ^ ( ' q q/2 q0 ^ to^Ti/2^ji4i9 :} Kp / )\/ 2A / ,^ 1 {lm d Im n \ Jm n + m d ) ( X0^m n 



4-7T 



= - q ) x , ,, r 



Im n 
(A20) 



where it is assumed that the incoming nucleons momentum is parallel to the z-axis, one gets the 
second matrix element (pqa\ t^j +c PGot^P |<J>): 



(pqa\4 +c PG t?P\<S>) = W fdp t^ R (p,p;E(q))G (p,q) 



1/2 1/2 t 



^ (« Pj E (00 I «o) C ao 6sS'$w' S M T ,-1/28 M T , -1/23V 
/3 <* A [ 1/2 T 1 

1/2 1/2 s 



y p ,2 dp^ ao ( P ') 



1/2 5 s' 



Xi(pqf3,pq f3 ) + 5 ss >X 2 (pq/3,p'q l3o) 



CA21) 



where |a >= |(Z / s')j"'(A'l/2)/'(//')^(*' 1 / 2 ) r ' = V 2 >, = 0,2, i' = 0, T' = 1 s' = / = 1 and 



C 



Q() 



— — — (lm d l'm n \ Jm n + m d ) ( \'0-m n 



1 m r 



(A22) 



APPENDIX B: PARTIAL WAVE PROJECTED TRANSITION MATRIX ELEMENTS: 

< pq\a' >< a'\T\$ > AND < pq\a' >< a'\PT\$ > 

Including spin (mj) and isospin (i/j) projections of the nucleons one gets for the kinematically 
complete breakup configuration specified by Jacobi momenta (p, q) the following contribution to 
the transition matrix element: 



( x pqm\rn2rnj,v\V2Vz\ / p' 2 dp' q' 2 dq' \p'q'a) (p'q'ot\ T |$) 



E E ^ M ~ mi " m2_m3 (w) E E 



J 71 £AL 



js/t S 



I S j 

A \ I 

L S J 
1 



(LSJ, M — mi — 771,2 — "7l + ?™2 + TTls) (^S^^' m 2 + 1^3,^1 

--s,m 2 ,m 3 \ (t^ T , v 2 + v$,v\) ( 2 2*' V ^ a l T l $ ) 



(Bl) 
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where the incoming state |3>) = |<fo, ^/x) \ifd, lrrid) is composed from the relative nucleon-deuteron 
motion with momentum % \\z, the deuteron wave function (p^, and M = fi + m d is a sum of spin 
projections of the incoming nucleon and deuteron. 

The corresponding contribution to the elastic scattering transition amplitude is: 

( $ 1 P Y1 J P' 2 *VW \p'q'a') (p'q'a'\ T|$) = 

a' 

J2 J2 ( 17 ° J > m 'di M - m 'd) ( A o^ J o, M - m' d - //, A Y X(hM _ m ^^(q' ) 

J"Ma>A /o ^ ' 

1 

/ q*dq' [ dxw {^) GaoAq ? q ' X) (^q'a'\T\^ (B2) 



where \a<i >= \(lol)l(Xo^)Io(Ho)J, (^)T = \ >, lo = 0,2, c^ ,2 is the deuteron wave function, and 



APPENDIX C: ($'| PGq 1 |$) AND MATRIX ELEMENTS WITH 3-DIMENSIONAL 
COULOMB T-MATRIX: (p<f m^l if P |$) AND ($'|PifP|$) 



Using 

P |p§*mim 2 m,3^i^ 2 ^3) = 
+ 



1 „ 3 ^ _ 1 ^ 
~2 P ~ 4 9 ' P_ 2 q,rn2m3mil/2iy3Ul / ' ' 



and 



1 1 \ 

'2 mjV 2 ? 

,M=m £ j— m2— 013 (p) (Lll.M^d-M) 



(p§mim 2 w, 3 ^i^2i / 3| = (p?mim2m 3 ^ii/2i / 3| VdXmdt = 00-m N -m t qoj = 



L=0,2 

--l,m 2 m 3 m d - Mj / --0, f 2 f 3 ) <5 mimjv <^ im( <5(<f - <f ) (C2) 



one gets the following expression for the contribution of the matrix element of the if P to the 
breakup transition amplitude 



(pqmm^P^) 



L 

(Lll, md — mi - m 2) rai + m 2 ) 
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(Lll,md — m 3 — mi,m 3 + mi) 



(C3) 



The corresponding contribution to the elastic scattering transition amplitude from Pt^P is 



($'|Pt*P|$) = 

£ £ I d^ L {\q+\q' \)Y^ mdi . 
mim 2 l,L 

(Lll\m,d> - mi - m 2 ,mi + m 2 , m' d )(-mi-m 2 |lmi + m 2 ) 
(Zll|m rf - mi -m 2 ,mi + m 2 , md)(-mi-m 2 |lmi + m 2 )<5. 



mi-m 2 (g + J? o)^i(l9+ 290l) y £,md-mi-ma(9+ 2«o) 



^(l 2 -<f+<f ol,||9 + ?o|,cos0 = 



|g 2 + go • go + |gj (go + g Q 

lk+9oll|g + ?ol 



...j. j. , ■■-^ / ~m N r,m N vmt N ,,m tN 



/l 11 
dq<p L (\q+ 7;? o\) Y L,m d ,-m 1 -m N (Q + ^ + ^\) Y L, 

"n l,L 

(Lll|m<f — mi — 771 7v, mi + mjv, m^)(-mi-mjv|lmi + mjv) 
(Lll|m<2 — m/v — mi, mjv + mi, m^)(-m7v/-mi|lmAf + mi)<5. 



2'--jv 2""- 1 - 1 J> ' ■•"- L '- m ijV' mt j V ' 

*/ i 1. 



o 1 w 1 jg +go-g Q+ 2 g" (gQ + g q) 3 2 

Ul^oU^'^ ^ +rom+ ^ ^ E -^ 



ere {p \t* \p ') = tf (\p \,\p '\,p ■ p ')• 

The matrix element < &'\PGq 1 | ( J> > is given by 



(C4) 



m p 4 



(- + cos9) + \E d \ 



LL'm,2 



go + ^go 



Y L' ,M'=m' d -rn N 



1 



:go + go 



- m2 (g + ^go ' ) (£'H, m' d - M'> / mjvm.m^ - M' 
)^L,M=m d -m 2 - miV (^go + g*o ' ) (Lll,M,m d - M) ^~l,m 2 m N m d - M^(C5) 



APPENDIX D: MATRIX ELEMENTS WITH 3-DIMENSIONAL COULOMB T- 
MATRIX: < pq mm\t*PGo\a' >< a'\T\$ > AND < &\Pt*PG \a' >< a'\T\$ > 

Inserting two complete set of states \pqmiVi > (mj and V{ are nucleons spin- and isospin- 
projections) and \pqa > one gets 



a; 



(pqmiV t \tKpG T\<5>} = J dp't*(p,p';E(q))J2 j dp" p" 2 dq" q" 2 (p ' qmm] P \p" q" c 

x G (p",q")(p"q"a\T\$) . (Dl) 
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The permutation matrix element {p'qmiVi\ P\p"q"a) is given by 

(p'qmiVi\ P \p"q"a) = £ / dp 4 dq 4 (p'qmiVi\P \p4q4mWi) {PiQifn'y^p' 1 'q" 'a 

E / dp 4 dq A 5{p' - vfi) 5 {p A + vf 2 ) 5 mim > 3 



+5{p' + 7Tl) 5 (p 4 - 1? 2 ) ^rn^rn'^rn^rn'^rn^rn^vrv'^v^Viv'^ 
I S j 

1 /I / 



E v 3 tL§ < 

LS 



/I ,1 , 








5(p»-p A )5(q"-q 4 )^ 


p//2 


q" 2 J 




2 m 2 -m 3 



>//;', + m'> ) ( sm' 2 + 771,3— 



SWj + 777 2 + m' 3 



L S J 



Th>[ + + f 3 \ (LM^Sm^ + m' 2 + m' 3 | Jm p + m^) 



(^4,54) 



where 



(D2) 



2* + 34 



^2 = g+2^4 



(D3) 



Performing in (|D2p integrations over and (/' leads to 

(pqmim 2 m 3 viV2V3\tcPGoT\®) = J dp'dp A dq A t^(p,p';E(q)) 



+6(p' + 7rl) (5 (p 4 - 1? 2 ) (Snur^^mam^msmi^^i/^^i/J 
/ s j 



Ev^< 



LS 




1 /I / 

777 9 — 777o 



Sm' 2 + 777 3 ^ ^ S777 2 + m'3- TO^ 



Sm'j + m 2 + m' 3 



1 ,1 , 

— Vn — Vn 

2 2 A 



Tf'i + u' 2 + f 3 ) (LMlStTIi + m'2 + 777 3 I JMj) 



^A Mi 94) G (f>4, 94) (P4, 94, a| T |$) 



(D4) 



For channels a ^ the only singularity is the Go singularity. When p 4 and q 4 are taken as an 
integration variables in (|D4|) Go becomes 

1 



G (P4,94) = G (p',q) 



(D5) 



£7 + te-±(p '2+3^2) 

and the contribution to the matrix element (pqmim2rn3UiU2i's\t^PGoT \^} coming from these 
channels is given by 

1 



(pg777 l7772777 3 ^i^2^3| t^PG T |$) | a ^ Qd = J p ' 2 dp ' 



E + ie 



p ' 2 3q 2 
m Am 



d(j) p 'd(cos 1 
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tf(p,p';E(q)) £ y/jft§< 

a^a d LS 



1 1 

-m 3 -mi 

1 1 

— ZA>— 

2 J 2 



57773 + mi ) { sm^ + m±-m2 



/ s j 

A | J J" {LM L Smi + 777,2 + "i3 I + m d ) 

L S J 
1 



^3 + v\ ) { tu 3 + v\-v 2 



Sm\ + 7772 + "73 
T^i + 7^2 + ^3 



IX 



2 P ~A q ' P ~2 Q 



3 1 

2^ - |9 I, \P ~ \' a 



T|$) 



, 1 1 
+ ( 2 mi 2 m2 



sm i — /??2 } { Smi + 777-2 — m 3 



Smi + 777-2 + 777-3 



1 1 

— 7/1 —7^2 

2 2 



*fl + ^2 ) ( + ^2^3 



Tui + + 



1^, 3. 
2^ 4 



/ t , 1-/ , 3„_, , ^ 1 



1^/ 3_ 
2 P + l ql 



P 



T|<&)] 



2" J \' T V '" ' 2 

where m p and are spin projections of the incoming proton and deuteron, and Ml = ra p + m d — 
mi — 7772 — m 3- Integration over p' can be performed numerically taking care of the pole in Go by 
e.g. the subtraction method. 

For a d channels one can decompose Go(j>4, (J4) (p494a t j|T|$) in (|D4|) in the following way 



(D6) 



^0(774,94) (p4q4a d \T\$) 



1 



4 m \Hmax 



m 



-Ma 2 

4m, Wr, 



^4 + ^4^(90"^)+^ 

! (p494ay|Tl$ 



(D7) 



Am^max Hi, ~P| + »C + ~P| ^(9o-«4)+* e 1^1 + ^4 

The first part in (|D7p can be integrated in the same way as the contribution from channels a ^ a d , 
resulting in expression (|D6jl with (pqa\T\<&) replaced by ^"^^^ ■ 

Calculating contribution from the second term in (1D7|) one takes in (|D3j) p ' and 7J4 as an 
integration variables. This gives the following result for that term 

1 



{pqmim2mzV\V2V2,\ t c PG T |<3?) 



I s j 



^i(9o-d)+* e 



dcf> q4 d(cos 9 q4 ) 



E v ^™ < 



A i J > (LM L Smi + 777 2 + 7773 j ^Tt p + 777 d ) 



L S J 



7 1 1 



S7773 + ?77l ) ( S7773 + 7771—7772 



1 1 

— 7/3 — 7/1 

2 2 



^3 + ^1 ) ( ^3 + ^1-^2 



Smi + 7772 + 7773 



1 



Ti/i + i/ 2 + ) 2« + 945 E(q)) 



*\{\-q-\q±\,q4,oLd\T\<&) 

~q - 2^4 , 94 — — — 2 — 

2 / +^9+594 
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II 1 

+ \2 mi 2 m2 



sm\ + rri2 



1 1 

— ZVi — V2 

2 2 




smi + 777-2—7773 



SYn-i + 771,2 + m 3 



Tvx + ^2 + ^3 ) if (p, -^g - 94! S(g)) 



1 



9+ 2941,94,^ \T |$) 



(D8) 



|^| + ^|9-+i94| 2 

Again one has a simple pole which can be treated using e.g. the subtraction method. 

For the exclusive breakup the calculation of the contributions (ID6I) and (|D8j) must be performed 
for each complete geometry in a coordinate system used when calculating < pqa\T\& > where the 
z- axis was taken to be parallel to the incoming proton momentum. 

For the elastic scattering an additional integration has to be done providing the following 
contribution to the elastic scattering amplitude from the < 3>'\Pt^PGo\a.' >< a'\T\$ > term: 



< &\Pt?PG \a' >< a'\T\<S> >= [- I dq 



(Lll,M,m' d -M,m' d ) 



mim 2 LM 

--l,m 1 m 2 m d - M 



• y L,Af(9+ 2% 



qo ,qmim2m N i,v 1 



1 



mim 3 LM 



1 

+ 2 mt N> 



1 

+ 2 



t?PG Q T\<f>) 



(Lll, M, m' d - M, m' d ) m 3 m im ' d - Al^j 



V2 
2 



-q + qo ', qmim N 'm3, v\ 



1 



1 

+ 2 



tfPG T|$) 



(D9) 
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FIG. 1: (color online) The convergence in the cut-off radius R of the pd elastic scattering cross section 4£ 
shown as a function of the cm. angle 8 c . m . at the incoming proton energy E l p ab = 13 MeV. These cross 
sections were calculated with the screened Coulomb force and the CD Bonn nucleon-nucleon potential 
restricted to the 1 5'o and 3 5i- 3 -Di partial waves. The screening radii are : R — 20 fm (black dotted line), 
R = 40 fm (green double-dashed -dotted line), R = 60 fm (blue long-dashed-dotted line), R = 80 fm (red 
dashed- double-dotted line), R = 100 fm (blue short-dashed line), R = 120 fm (red long-dashed line), 
R = 140 fm (black solid line). The R = 40-140 fm lines are practically overlapping. The black dashed- 
dotted line is the corresponding nd elastic scattering cross section. The pluses are E l £ h — 12 MeV pd elastic 



scattering cross section data of Ref. 
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FIG. 2: (color online) The convergence in the cut-off radius R of the pd elastic scattering cross section 4£ 
at the incoming proton energy E l p ab = 13 MeV, shown as the ratio ^ R /^ R 14 °- For the description of the 
lines see Fig. Q~] 
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FIG. 3: The contributions of different terms to the pd elastic scattering cross section at the incoming pro- 
ton energy E l p ab = 13 MeV calculated with the screening radius R = 100 fm. The dotted and dashed lines are 
contributions of the ($'|PT|$) and ($'|PGq |$) terms, respectively. The double-dashed-dotted line is the 
contribution of the ($'|P^P|$) term coming with the 3-dimensional screened Coulomb t-matrix t^. The 
dashed-double-dotted and dashed-dotted lines are contributions of the ($'|PtfP|$) and (&\PtfPG T\§) 
terms, respectively, which are calculated with the partial- wave decomposed screened Coulomb t-matrix. The 
solid line encompasses all terms. In this feasibility study the 3-dimensional tf is replaced by Vj 1 . 
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FIG. 4: (color online) The independence of the real (left column) and imaginary (right column) parts of 
the partial wave contribution (®' md , mpl \P(T - tfP - tf PG T) | $ m<jTOp ^ to the E l p ab = 13 MeV pd elastic 
scattering transition amplitude on the cut-off radius R. The different lines are: R = 20 fm - dotted, R = 40 fm 
- short-dashed, R = 60 fm - long-dashed, R = 80 fm - short-dashed-dotted, R — 100 fm - long-dashed-dotted, 
R = 120 fm - double-dotted-dashed, R — 140 fm - solid. All the lines are practically overlapping. The 
incoming and outgoing deuteron and proton spin projections are for a): = = — 1 and m p = m p i = — \ 
and for b): m,d = — 1, mji = +1, and m v = m p > = +i. 
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